Spectral/hp least-squares finite element formulation for the Navier–Stokes equations

We consider the application of least-squares finite element models combined with spectral/hp methods for the numerical solution of viscous flow problems. The paper presents the formulation, validation, and application of a spectral/hp algorithm to the numerical solution of the Navier–Stokes equations governing two- and three-dimensional stationary incompressible and low-speed compressible flows. The Navier–Stokes equations are expressed as an equivalent set of first-order equations by introducing vorticity or velocity gradients as additional independent variables and the least-squares method is used to develop the finite element model. High-order element expansions are used to construct the discrete model. The discrete model thus obtained is linearized by Newton’s method, resulting in a linear system of equations with a symmetric positive definite coefficient matrix that is solved in a fully coupled manner by a preconditioned conjugate gradient method. Spectral convergence of the L₂ least-squares functional and L₂ error norms is verified using smooth solutions to the two-dimensional stationary Poisson and incompressible Navier–Stokes equations. Numerical results for flow over a backward-facing step, steady flow past a circular cylinder, three-dimensional lid-driven cavity flow, and compressible buoyant flow inside a square enclosure are presented to demonstrate the predictive capability and robustness of the proposed formulation.     

A pseudo‐viscoelastic approach for transient polymeric flow exiting a channel

The influence of elasticity of a fluid exiting a channel is examined on transient coating downstream. A hybrid spectral/boundary element approach is proposed to solve the problem. The flow inside the channel is assumed to be fully developed. A viscoelastic instability of one‐dimensional plane Couette flow is first determined for a large class of Oldroyd fluids with added viscosity, which typically represent polymer solutions composed of a Newtonian solvent and a polymeric solute. The Johnson–Segalman equation is used as the constitutive model. The velocity profile inside the channel is taken as the exit profile for the emerging free‐surface flow. The flow is assumed to be Newtonian as it emerges from the channel. An estimate of the magnitude of the rate‐of‐strain tensor components in the free‐surface region reveals that they are generally smaller than the shear rate inside the channel. The evolution of the flow front is simulated using the boundary element method. For the channel flow, the problem is reduced to a nonlinear dynamical system using the Galerkin projection method. Stability analysis indicates that the channel velocity may be linear or non‐linear depending on the range of the Weissenberg number. The evolution of the coating flow at the exit is examined for steady as well as transient (monotonic and oscillatory) channel flow. It is found that adverse flow can exist as a result of fluid elasticity, which can hinder the process of blade coating. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of the unsteady behaviour of cavitating flows

A 2D numerical model is proposed to simulate unsteady cavitating flows. The Reynolds‐averaged Navier–Stokes equations are solved for the mixture of liquid and vapour, which is considered as a single fluid with variable density. The vapourization and condensation processes are controlled by a barotropic state law that relates the fluid density to the pressure variations. The numerical resolution is a pressure‐correction method derived from the SIMPLE algorithm, with a finite volume discretization. The standard scheme is slightly modified to take into account the cavitation phenomenon. That numerical model is used to calculate unsteady cavitating flows in two Venturi type sections. The choice of the turbulence model is discussed, and the standard RNG k–εmodel is found to lead to non‐physical stable cavities. A modified k–εmodel is proposed to improve the simulation. The influence of numerical and physical parameters is presented, and the numerical results are compared to previous experimental observations and measurements. The proposed model seems to describe the unsteady cavitation behaviour in 2D geometries well. Copyright © 2003 John Wiley & Sons, Ltd.     

A computational method for the hydrodynamics of fractured-porous media

A numerical method based on the boundary-fitted finite difference method (BFDM) is presented in this paper. The boundaries are external (the boundary of the physical domain) and internal (which corresponds to the fracture network). The difference between this approach and the usual one lies in the inclusion of discrete fractures in the volume that represents the porous medium. The numerical model has been used in the prediction of the flow pattern in several internationally recognized verification cases and applied to the solution of hypothetical problems of interest to us in the field of nuclear waste repository modelling. The results obtained show that the numerical approach considered gives accurate and reliable predictions of the hydrodynamics of fractured-porous media, thus justifying its use for the above-mentioned studies.     

考虑行车状态的一维元胞自动机交通流模型

在Nagel Schrekenberg单车道元胞自动机交通流模型（简称NS模型）的基础上，考虑车辆之间的相对运动薛郁等提出了一种改进的单车道元胞自动机交通流模型（简称改进的NS模型）.通过两种情况列出了改进的NS模型存在不尽周严的地方，随之在新模型中引入了行车状态 变量和反馈规则，从而控制车辆出现倒车和刹车过急等现象.通过计算机对新模型进行模拟 ，发现减速概率和车流密度对车流状态的演化影响很大，当减速概率高（如道路条件差）时 ，即使车流密度低，车流也会出现局部堵塞状态；而当减速概率一定时，随着车流密度增加 ，车流的运动相与堵塞相发生了全局性的交替出现，此时类似于波的波峰和波谷的传播.与 改进的NS模型相比较，新模型模拟的车流量较高，说明新模型减少了车流的总体停滞状态. 关键词： 交通流 元胞自动机 行车状态 反馈规则     

Optical resonator for high-power transverse flow CO2 lasers

A new design of the U-type resonator is described. In this way, a laser beam with symmetrical intensity profile (regarding to a symmetry plane) can be extracted from an active medium that exhibits gain asymmetry along one of the transverse directions. The whole area of the active medium cross-section can be used, and consequently the laser efficiency will be increased. This resonator structure was applied for efficiency power extraction (as a low order TEM modes laser beam) from a DC excited transverse flow CO₂ laser with cylindrical geometry. Although the cross-section area of the discharge was entirely used (including the cathode fall region), a symmetrical intensity profile of the laser beam (regarding to the two orthogonal symmetry planes) was obtained in the near field as well as in the far field; the gain asymmetry along the flow direction was compensated by the gas circulation fluidodynamical circuit with two counterflowing discharge channels. A double-U optical resonator was introduced in order to provide a laser beam with axial symmetry.For the practical construction of these two types of optical resonators we have developed two new types of 90° deflection elements: the first one, which does not reverse the image (and which has the properties of the pentaprism), and the second one, which rotates the image with 90° angle. Both elements exhibit good focusability if they are equipped with two concave mirrors.     

双侧进气突扩燃烧室内液雾燃烧的数值模拟

本文采用全双流体模型对部分切向进气管式双侧进气突扩燃烧室内的液雾燃烧过程进行了数值模拟，将液雾群按初始尺寸分组，数值模拟结果表明不同尺寸液雾的蒸发及燃烧过程显著不同，小尺寸的液雾是先蒸发后燃烧，而较大尺寸的液雾则是边蒸发边燃烧．模拟结果还表明，部分切向进气方式有利于在燃烧室的头部产生大面积的中心回流区，该回流区的存在有利于加强混合及火焰稳定．     

Chaos,transport and mesh convergence for fluid mixing

Chaotic mixing of distinct fluids produces a convoluted structure to the interface separating these fluids. For miscible fluids （as considered here）, this interface is defined as a 50% mass concentration isosurface. For shock wave induced （Richtmyer-Meshkov） instabilities, we find the interface to be increasingly complex as the computational mesh is refined. This interfacial chaos is cut off by viscosity, or by the computational mesh if the Kolmogorov scale is small relative to the mesh. In a regime of converged interface statistics, we then examine mixing, i.e. concentration statistics, regularized by mass diffusion. For Schmidt numbers significantly larger than unity, typical of a liquid or dense plasma, additional mesh refinement is normally needed to overcome numerical mass diffusion and to achieve a converged solution of the mixing problem. However, with the benefit of front tracking and with an algorithm that allows limited interface diffusion, we can assure convergence uniformly in the Schmidt number. We show that different solutions result from variation of the Schmidt number. We propose subgrid viscosity and mass diffusion parameterizations which might allow converged solutions at realistic grid levels.     

G-spaces and Fixed Point Sets

For a compact connected group G, in a compact free G-space M every closed subset is the fixed point set of an autohomeomorphism if the weight of M is not bigger than that of G. In a compact group with nonsingleton components, every closed subset is the fixed point set of an autohomeomorphism if and only if the weight of the whole group is not bigger than that of the component of the identity.     

RHIC能区Au+Au碰撞中集合径向流的研究

利用相对论量子分子动力学模型RQMD,对RHIC能区s=200A GeV Au+Au碰撞的集体膨胀效应进行了研究,对散射粒子的横质量谱进行了分析.研究表明,在RHIC能区的重离子反应中存在有强的集合径向流.对单粒子谱的拟合结果给出Au+Au的源冻结温度为160MeV,平均径向流速度为0.6c.